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In  this  paper,  we  present  an  Age-Structured  hepatitis  B  model.  This 
epidemic  model  is  investigated  for  different  classes  of  infectious  diseases 
that  can  be  transmitted  through  an  effective  contact  with  infective 
individuals,  who  are  contagious.  The  Graph-Theoretic  Method  for  the 
Basic  Reproduction  Number  was  obtained.  In  addition,  the  numerical 
simulation  is  used  to  virtually  verify  the  model  predictions.  The  result 
suggest  that  the  endemic  nature  of  the  model  is  approaching  equilibrium 
with  increase  immunization  program  and  other  control  measures  put  in 
place. 


1.  Introduction 

In  the  last  two  decades,  mathematical  models  have  frequently  been  used  to  study  the  transmission 
dynamics  of  HBV  in  several  regions.  Anderson  and  May  [1]  used  a  deterministic,  compartmental 
mathematical  model  to  illustrate  the  influences  of  carriers  on  the  transmission  of  Hepatitis  B  virus 
(HBV).  Anderson  et  al.  [2]  and  Williams  et  al.  [3]  described  models  of  the  sexual  transmission  of 
HBV,  which  include  heterogeneous  mixing  age  and  sexual  activity.  Edmunds  et  al.  [4]  illustrated 
the  relation  between  the  age  at  infection  with  HBV  and  the  development  of  the  carrier  state.  In 
2007,  Stanca  et  al.  [5]  modeled  the  mechanisms  of  acute  hepatitis  B  virus  infection,  they 
discovered  that  a  cell-mediated  immune  response  plays  a  significant  role  in  controlling  the  virus. 
Thomley  et  al.  [6]  applied  the  model  of  Medley  et  al.  [7]  to  predict  chronic  hepatitis  B  infection  in 
New  Zealand.  Hepatitis  B  is  a  severe  liver  disease  caused  by  the  hepatitis  B  virus  (HBV),  Edmund 
et  al.  [4].  It  is  a  primary  global  health  problem  and  the  most  severe  type  of  viral  hepatitis. 
Fonnerly  known  as  ’’serum  hepatitis,”  the  infection  has  caused  epidemics  in  parts  of  Africa  and 
Asia;  moreover,  it  is  endemic  in  China,  Williams,  [8]. 

2.  Mathematical  Model  Formulation 

We  propose  a  mathematical  model  to  understand  the  transmission  dynamics  and  prevalence  of 
HBV.  The  model  is  constructed  based  on  the  characteristics  of  HBV  transmission  and  the  model 
of  Sirajo  et  al.  [9]. The  model  comprises  of  nine  compartments,  here  we  see  the  need  of  treatment 
of  carrier  above  15  years  of  age,  because  of  the  need  to  procure  a  therapeutic  treatment  for  the 


158 


A.K.  Dotia  et  al.  /  Journal  of  Science  and  Technology  Research 
1(1)  2019  pp.  158-165 

infected  individuals  in  the  population.  Susceptible  individual  below  15  years  of  age  ( S(J  (7) ), 

Susceptible  individual  15  years  and  above  (SF  (7) ),  Vaccinated  V,  Infectious  individual  below  15 

years  of  age  (/^(r)),  Infectious  individual  15  years  and  above  ( IF(t )),  Chronically  infected 

individual  below  15  years  of  age  ( Cu  (f)  ),  Chronically  infected  individual  above  15  years  of  age  ( 

CF  (t)  ),  Treatment  of  chronically  infected  individual  above  15  years  of  age  and  above  (  Tp  (f) ), 

Recovered  (R(t)),  and  Vaccinated  (V(t)).  The  following  assumptions  were  made:  There  exist 
disease  induced  deaths  due  to  chronic  (via  single  and  dually  HBV/HIV  co -infections).  Treatment 
should  be  for  a  long  time  and  uninterrupted,  as  virological  relapses  after  discontinuation  of 
treatment  are  frequent.  Wanning  of  HBV  vaccination  takes  25  years,  and  it  has  not  been 
established  to  be  life-long. 


,T,V(1— 

1 


i  1  1 


Figure  1:  Schematic  representation  of  interactions  of  HBV  transmission 
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—  rccF  _  (/^o  Ri +/t)tf 


dt 
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(1) 

(2) 
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(4) 

(5) 

(6) 

(7) 

(8) 
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Where, 

.  n  =  epTb ,  0)f  =  £pTf  ,coc  =  SCTC  ,nA  =  7TV  (1  -  coh ) ,  nc  =  KV(j)  (1  -  0Jh  ),PA=pc(\-cah), 

Pc  =  pcp{\-(oc),rAu  ~  ^ A  (l  dn)  ^  Y  AF  =°a{]-  df )  ,  Z0  =  p0  +  crA,Zi  =  p0  +  cr  +  cou ,  Z2  =  p()  +  cof , 
Z3  —  p0  +  y/ ,  Z4  —  p0  +  p1+  o' A ,  Z5  —  p0  +  <7 c ,  Z6  —  p0  +  P2  +Yc>  ^7  —  Mo  Mi  Yt  > 


Table  1:  The  list  of  parameters  /  notations 


Parameter 

Interpretation 

Value 

Reference 

M 

birth  rate 

0.0121 

[11] 

Mo 

Natural  mortality  rate 

0.00693 

[11] 

Mi 

HBV  related  mortality  rate  by  I U  and  I F 

0.007 

[10] 

Mi 

HBV  related  mortality  rate  by  C  F  and  T  F 

0.00131 

[9] 

¥ 

Rate  of  waning  of  vaccine-induced  immunity 

0.04 

[9] 

q 

Average  probability  an  individual  fail  to  clear  an  acute  infection 
and  develops  to  carrier  state 

0.885 

[9] 

n 

birth  rate 

0.036 

[9] 

c 

Average  number  of  sexual  partner 

3.233 

[9] 

p 

HBV-Sexual  transmission  risk  rate 
and  pc  is  the  effective  contact  rate 

0.6 

[9] 

d 

Modification  Parameter  that  suggest  reduce 
sexual  transmission  rate  by  chronic  individual 

0.667 

[9] 

Condom  efficacy 

0.8 

[9] 

£P 

Vaccine  efficacy 

0.9 

[9] 

V 

Proportion  of  perinatal  infected  HBV  positive  birth 

0.724 

[91 

<f> 

Modification  Parameter  that  suggest  reduction 
in  HBV-positive  birth  by  chronic  individual 

0.159 

[9] 
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Rate  of  moving  from  S  U  to  S  F 

0.0667 

[9] 

® A 

Rate  of  moving  from  acute  to  chronic  infection 

2.667 

[9] 

ac 

Rate  of  moving  from  C  U  to  C  F 

0.069 

[9] 

Yt 

Rate  of  moving  from  T  F  to  R 

Hypothetical 

Yc 

Rate  of  moving  from  C  F  to  T  F 

0.015 

[9] 

Clu 

Proportion  of  I  U  which  progress  to  C  U 

0.885 

[12] 

qf 

Proportion  of  I F  which  progress  to  C  F 

0.1 

[9] 

TC ’  rfc’  T«»  Tf 

Condom,  Vaccine(birth),  Vaccine 
(below  15),  Vaccinefabove  15),  compliance 

0.1 

[9] 

The  total  population  N(t)  can  be  obtained  from 

N(t)  =  Su(t)  +  Sp(t)  +  V(t)  +  Iu  (t)  +  IF  (t)  +  Cv(t)  +  CF  (t)  +  TF(t)  +  R(t) 

^-  =  kN-  jU0N  -jU2(Cf+Tf) 

N(t)  =  N(  0)e{^o)'  (10) 


7T 


in  the  absence  of  disease  N(t)  — >  — ,  Moreover,  under  the  dynamics  described  by  the  above 

Mo 

systems  of  equations,  the  region 

\x  =  (Su,SF,V,Iv,Ip,Cu,Cr,TF,N)e*R9+\Su>0,  SF  >  0,  V  >  0,  Iv  >  0,  IF>0,Cv>0, 


cf>o,tf>o,n<^- 

Mo 


(11) 

is  positively  invariant.  Hence  the  system  is  both  mathematically  and  epidemiologically  well- 
posed.  Therefore,  for  initial  starting  point  x  e  S.K9 ,  the  trajectory  lies  in  Q.  Thus  we  restrict  our 
analysis  to  the  region  Q.  (where  the  models  make  biological  sense) 


2.1  Positivity  of  Solutions 

Lemma  1:  All  the  solution  of  the  Equations  1-9  are  positive  for  all  time  t  >  0  provided  the  initial 
condition  are  positive. 

Proof:  Let  {($„  (0) ,  SF  (0) ,  V  (0) ,  Iv  (0) ,  I F  (0) ,  Cv  (0) ,  CF  (0) ,  Tp  (0) ,  R  (0))  >  o}  e  9f+ 

C^-  =  7rN(l-cob)-7rv(lF+(/>CF)(l- 


')Sl j  Ov  +  cou  +  //(l )  SL, 

{ 12 ) 

)SV  ~(as+co„  +  Mo ) Su 

(13) 

This  implies  S[;  (f)  >  —jtv (lF  +  (j)CF  ) (l - coh ) Su  -(oy  +  cou  +  //„ ) Su  ,  integrating 
we  have 

Su(t)>Su{0)e  3  >0 


(14) 
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similarly  it  can  be  shown  that  Sv  >  0,  SF  >  0,  V  >  0, /,,  >  0,  IF  >  0,  Cv  >0,CF>0,TF>0,R>  0,  for 

all  time  t  >  0.  Hence  all  solutions  of  the  HBV  model  remain  positive  for  all  non-negative  initial 
conditions. 


3.  Basic  reproduction  number  (Graphic-Theroretic  Method) 

A  graph-theoretic  method  for  calculating  R0  was  discussed  exclusively  in  Tomas  et  al.[13],  [14]. 
From  the  definition  ofi?,,  =  p(FV  '),  they  are  able  to  derive  a  series  of  rules  for  reducing  the 
digraph  associated  with  FT-1  -V  to  a  digraph  with  zero  weight,  from  which  /.  =  R0  is  given.  The 
rules  are  as  follows: 

Rule  1.  To  reduce  the  loop  -  an  <  0  to  -  1  at  node  i,  every  arc  entering  i  has  weight  divided  by  an  . 
Rule  2.  For  a  trivial  node  i  on  a  path  j  — >  i  — >  k,  the  two  arcs  are  replaced  by  j  — >  k  with  weight 
equal  to  the  product. 

Using  Graph  Theoretic  Method 
Taking  the  infective  classes  of  the  model; 


dt 


7tv  [i  F  +  (ffC  F^(\  (ob)Su  (aA+pQ)IL 


dIF  _  pc(lF+TjCF)(\-ac) 


dt 

dCv 

dt 

dCF 

dt 


N 


+  R1+°a)If 


aACtu  Ijj  (  °C  A)  )  C(J 

—  °CCU  _  ( A)  J"  A?  Yc  )  C ! 


We  have  four  infected  classes  in  this  model,  Infectious  individual  below  15  years  of  age  ( Iv  (t) ), 
Infectious  individual  15  years  and  above  (IF  (r)),  Chronically  infected  individual  below  15  years 
of  age  ( Cv  (r)  ),  Chronically  infected  individual  15  years  and  above  (  Ch  ( t )  ),  hence  our  m  =  4. 

(15) 


0 

Z0 /u  KJfSU  ncCFSU 

F  = 

{$Jf  +  $cCf)sf 

,  v= 

ZJu 

0 

Z5CU  ~  °AquU 

0 

Z(f  l  ~  °Aq/F  ~  CCCU 

o 

o 

o 

o  1 

r  z« 

F  = 


0  PAsf 


0 

0 


0 

0 

0 


PcsF 

0 

0 


(16) 


V  = 


0 

~aAClu 

0 


71 A$U 


0 


-°AClf 


0 

0 

-<?r 


XcSu 

0 

0 


(17) 


For  the  disease  free  equilibrium  point  of  the  system  of  equations,  which  has  the  coordinates: 

s*  _nN(\-a>b)  s*  = _ Z4(-qu7tcCTACTc+ZbZ5ZA) _ 

Z\  n AquPco a<jc  —  7VcqupA(TA(Tc  +  Z4ZAq  fPcoA  +  fiAZ(ZZZ4 
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FA-1- 


R0- A 


Z0 

0  /3aSfA1  -Z4 

0 

0  oAqf 


0  ncSv 

0  pcsFx _1 

-Z5  0 

°c  ~Ze 


(18) 


Creating  A.  Digraph  FX  1  —  v 


f  PcG.Sf^F^ 

nc^u  j 

\zt{zi-pJsFri) 

k  Z6  z6(z4-/sa^1)> 

equating  the  left  hand-side  to  zero,  we  have; 
Pc^a^f  {^a^u^iFc  +  Z5Z0qf  )  a coAq  un  CS  v  _ 

Z(A  ( — fiAS  F  +  Z4  )  Z5Z0  Z6Z5Z0 

5F  f(^/d„  {~KaPc  +  ^cPa  )  -  PcZ,ZSf  )  ^  -PaZ6Z5Zo) 

Z4  {^C<JC(7AClu Su  ~  Z6Z5Z0  ) 


(19) 

(20) 
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4.  Results  and  Discussion 

We  simulated  numerically  with  the  aid  of  Maple  18  software  [15]  to  check  and  determine  the 
effect  and  behavior  of  the  parameters  of  the  model.  This  illustrate  the  asymptotic  stability  of 
Disease  free  equilibrium  (DFE)  of  the  model.  The  model  described  by  Equations  (1-9)  exhibit  a 
rich  dynamic.  We  observed  that  varying  these  control  parameters;  ib,  xu,  tc,  if,  when  there  is  no 
control  (effective  immunization’s  and  condom  usage  equal  zero)  our  RQ  =  2.9479,  at  low  rate  (10 

percent)  RQ  =  1.6145,  at  moderate  rate  (50  percent)  Rtl  =  0.3823,  at  high  rate  (90  percent)  Rtl  = 

0.0987.  However,  the  table  below  also  show  the  gradual  decrease  in  the  reproduction  number 
when  different  control  measures  put  in  place  are  increased.  The  result  established  here  is  more 
effective  in  bringing  the  basic  reproduction  number  below  unity,  when  compared  with  the  results 
from  literature  owing  to  the  treatment  class  adopted. 


Table  2:  Basic  Reproduction  Number  ( R0 )  obtained  for  HB  V 


Parameter 

RTb 

RTo 

Rr 

'Ht 

Rf 

JV0 

Tb,  Tu,  Tc, 

Tf 

DFE 

Remark 

DFE 

Remark 

DFE 

Remark 

DFE 

Remark 

0.1 

0.4664 

Stable 

0.3767 

Stable 

0.6115 

Stable 

1.2621 

Unstable 

0.2 

0.4454 

Stable 

0.3795 

Stable 

0.5752 

Stable 

0.8847 

Stable 

0.3 

0.4244 

Stable 

0.3809 

Stable 

0.5389 

Stable 

0.6811 

Stable 

0.4 

0.4033 

Stable 

0.3817 

Stable 

0.5026 

Stable 

0.5536 

Stable 

0.5 

0.3823 

Stable 

0.3822 

Stable 

0.4663 

Stable 

0.4664 

Stable 

0.6 

0.3612 

Stable 

0.3827 

Stable 

0.4301 

Stable 

0.4029 

Stable 

0.7 

0.3400 

Stable 

0.3830 

Stable 

0.3938 

Stable 

0.3546 

Stable 

0.8 

0.3188 

Stable 

0.3832 

Stable 

0.3576 

Stable 

0.3166 

Stable 

0.9 

0.2976 

Stable 

0.3834 

Stable 

0.3213 

Stable 

0.2860 

Stable 

1.0 

0.2763 

Stable 

0.3835 

Stable 

0.2850 

Stable 

0.2608 

Stable 

5.  Conclusion 

In  conclusion,  the  results  suggests  that  the  more  we  vaccinate  and  provide  adequate  treatment  to 
the  population  and  also  enlighten  them  on  its  menace,  this  will  bring  the  threshold  value  below 
unity,  thus  eradicating  the  disease.  Moreover,  other  highly  sensitive  parameters  are  ca  ,  cc  ,  as  ,q, 
yc  and  ytf.  For  the  parameters  a  a  and  qf  are  the  movement  from  an  acute  stage  to  a  chronic  stage 
and  should  be  of  high  sensitivity.  The  acute  stages  have  short  duration  (six  months)  while  chronic 
stages  have  the  long  duration  of  years  and  in  most  cases  they  are  life-long.  For  q,  it  is  the 
modification  parameter  that  suggests  the  reduced  sexual  transmission  rate  by  chronic  individuals, 
and  thus,  it  profoundly  affects  the  transmission  dynamics  of  the  disease. 
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